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The information about the current-carrying ability of a type-II superconductor can be obtained 
by studying the flux relaxation and vortex penetration phenomena in the superconductor. In early 
studies, the infinite series models of the flux relaxation and vortex penetration phenomena were 
constructed at a vanishing current density and vanishing internal field, respectively. However, this 
is not the only possibility. Here it is shown that one can reconstruct the theoretical models at the 
critical points. The new polynomial model of the flux relaxation (vortex penetration) phenomenon 
was constructed by expanding the vortex activation energy as an infinite series of the current density 
(internal field) about the critical current density (equilibrium internal field). The unification of the 
polynomial models was proposed. The inverse model of the flux relaxation (vortex penetration) 
phenomenon was also constructed by expanding the vortex activation energy as an infinite series 
of the inverse current density (inverse internal field) about the critical current density (equilibrium 
internal field). 

PACS numbers: 74.25.Op, 74.25.Uv, 74.25.Wx 



I. INTRODUCTION 

In the study of a flux relaxation and vortex penetra- 
tion phenomenon, an important step is to find out an 
expression of the vortex activation energy. The flux re- 
laxation and vortex penetration can be regarded as the 
processes of vortex hopping between the adjacent pin- 
ning centers fir— Therefore, they can be described by the 
Arrhenius relation 2 ^^, which shows that the vortex hop- 
ping rate has a strong dependence on the temperature 
and vortex activation energy. Because the vortex hop- 
ping rate can be changed by a current density (or internal 
field), the activation energy is a function of the current 
density (or internal field). By proposing a detailed activa- 
tion energy, one can obtain the corresponding time evolu- 
tion equation of the current density (or internal field) . A 
number of activation energies were proposed on the basis 
of different physical considerations. In the flux relaxation 
process, we have linear activation energy^, quadratic ac- 
tivation energy, series activation energy^, inverse-power 
activation energy^ and logarithmic activation energy&£. 
In the vortex penetration process, we have linear activa- 
tion energy, quadratic activation energy, and series acti- 
vation energy?^ 

One can see that the flux relaxation and vortex pene- 
tration phenomena are currently described with multiple 
models. The present work is motivated by the following 
facts: 

First, it is now understood that the early devel- 
oped infinite series models of flux relaxation and vor- 
tex penetration^— are constructed at zero points (vanish- 
ing current density and vanishing internal field). From 
mathematics, we know that one needs to choose a "cen- 
ter point" when doing a Taylor series expansion of a 
function. This means that it is possible to reformu- 
late the theoretical models by choosing the new center 
points. The critical current density and equilibrium in- 



ternal field are important parameters in the application 
of a superconductor,^ Therefore, we are interested in 
reconstructing the polynomial models at these critical 
points. 

Second, the inverse-power activation energy^ and log- 
arithmic activation energy^ are functions of inverse cur- 
rent density. They give good description to the flux re- 
laxation process over the random pinning in the high-T c 
superconductors. This led us to consider whether the "in- 
verse current dependence" implies a more general phys- 
ical law that includes the inverse-power activation en- 
ergy and logarithmic activation energy as special cases. 
It could be meaningful to extract out this common prop- 
erty and describe it with a unified mathematical equation 
which can be applied to a wide class of vortex systems. 

Finally, the activation energy of a vortex penetration 
process is currently expressed as a function of the in- 
ternal field. 2,3 But the activation energy with an inverse 
internal field dependence is unavailable. Similar to that 
in the flux relaxation process, therefore, it is desirable 
to propose an inverse model for the vortex penetration 
process. The inverse models must be constructed at the 
critical points because they diverge at vanishing current 
density (or vanishing internal field). 

In this article, we expanded the vortex activation ener- 
gies at the critical points and developed the new polyno- 
mial models and inverse models of vortex dynamics. Sur- 
face barrier— is included in the activation energy, but ge- 
ometrical barrier— li 3 - is ignored because it depends on the 
shape of sample. The information about the anisotropy 
of the superconductors is included in the activation en- 
ergy, which is related to the coherence length and pene- 
tration depth. 
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II. CLASSIFICATION OF VORTEX 
ACTIVATION ENERGY 

A number of vortex activation energies were proposed. 
Although these activation energies were constructed 
on the basis of different physical considerations, their 
mathematical expressions show common features. By 
carefully checking the existing activation energies, we 
find that they fall into two classes: 

Class A: Polynomial Activation Energy 

In the flux relaxation process, we have: 

1. Series activation energy^ U r (j) = Uo — Ya=i a tf l , 

2. Quadratic activation energy^ U r (J) = Ua—aij—a2j, 

3. Linear activation energyi^ U r (j) = Ua — a\j. 
These activation energies are decreasing functions of 

the current density j and are well-defined functions for 
all j G [0, j c ], where j c is the critical current density. 
In the vortex penetration process, we have: 

1. Series activation energy^ U P (B) =Uq + J2iZi c iB l , 

2. Quadratic activation energy U P (B) — U0+C1B+C2B 2 , 

3. Linear activation energy U P {B) = Uq + c\B. 

These activation energies are increasing functions of 
the internal field B and are well-defined functions for 
all B G [0,B e ], where B e is the equilibrium internal 
field (maximum internal field) of the vortex penetration 
process. 

Class B: Inverse Activation Energy 

In the flux relaxation process, we have: 

1. Inverse-power activation energy^ 7 U r (j) = Uo(j c /j)^, 

2. Logarithmic activation energy 8 '- U r (j) = Uoln(j c / j). 
These activation energies are functions of the inverse 

current density 1/j and diverge at j = 0. In the vor- 
tex penetration process, the inverse activation energy is 
currently unavailable. 

In Class A, the infinite series activation energies are 
the generalizations of other activation energies. From 
mathematics, we know that the Taylor series of a function 
f(x) about a point x$ can be written as 

00 

f(x) = f[x + (x- x Q )} = ^ a i( x - x oY, (!) 

1=0 

where a; = f^(xa)/ll. This shows that, to express a 
function as a Taylor series, one has to choose a center 
point Xo. By carefully checking the infinite series activa- 
tion energies in Class A, we see that they are the Taylor 
series expansions about the point xo = 0, more specif- 
ically, the expansion of U r about the point j ' = and 
Up about the point B = 0. These series are known as 
Maclaurin series. However, there is no reason why we 
limit ourselves to the zero point. We should consider the 
possibility of constructing the new polynomial models of 
the flux relaxation and vortex penetration by choosing 
other center points. The critical points are good candi- 
dates. 



On the other hand, the activation energies in Class B 
have inverse current dependence. These activation ener- 
gies can be regarded as the special models constructed at 
the critical current density j c . Therefore, we should gen- 
eralize this "inverse current dependence" to an infinite 
series and apply this method to the vortex penetration 
phenomenon as well. 



III. POLYNOMIAL MODEL 

In this section, we first constructed a new polynomial 
model of flux relaxation by expanding the activation en- 
ergy U r about the critical current density j c . Next, we 
constructed a new polynomial model of vortex penetra- 
tion by expanding the activation energy U p about the 
equilibrium internal field B e . 



A. Polynomial model of flux relaxation 
phenomenon 

1. Current dependence of the activation energy in a flux 
relaxation process 

In the flux relaxation process, the activation energy 
U r is a decreasing function of the current density j ~£ To 
construct the polynomial model at the critical current 
density j c , let us define a normalized current density a = 
j/ jc- We have a G [0, 1] because j G [0, j c \. 

Referring to Eq.((T]) and rewriting the activation energy 
as U r [l + (a — 1)], we see that U r needs to be expanded 
as a decreasing function of (a — 1) about the point a = 1. 
Let us explicitly write out the series expression of U r as 

00 

tf r (a) = -$>,(«-!)', (2) 
1=1 

where U r (l) — (U r is zero at j c ), ai — —U' r (l), 0,2 = 
-U?(l)/2\,---,a l = -U r l) (l)/l\. 

With Eq. @ , one can directly go to Subsection IIII A 21 
and derive the time dependence of the current density. 
Before we do that, let us discuss the pinning potential 
and the vortex activation energy at a vanishing current 
density. 

In a flux relaxation process, the external field is usu- 
ally zero. At a vanishing current density, a vortex inside 
the superconductor is subjected to a pinning force. If 
the vortex is close to the surface, then it is further sub- 
jected to an attractive surface imaging force. This means 
that the vortex activation energy at the vanishing current 
density, U r o, should includes the pinning potential U c and 
the reduction to the activation energy caused by the sur- 
face imaging force Ui m = / AirX) 2 K${2x / X) . Thus, we 
have 
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The coefficients bi in Eq. fTO]) are 1 



(3) 



On the other hand, U r o can be obtained by putting 
a = in Eq.fl2Jl, that is, 



Oj. 



(4) 



Combing Eq.© and Eq.flU), we have 



;=i 



2x 



The coefficients a/ can be obtained from experiments 
(see Eq.(|12p). Doing a measurement at a position away 
from the surface where Ui m — / Att\) 2 K (2x / \) = 0, 
one can calculate the pinning potential U c using Eq.©. 



h =37 E (-D s+i+u+ - 



l l 
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l{l+l)---(l-l + s + t + u + ---) 



s\t\u\ ■ ■ 

t 
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(11) 



where s + 2t + 3u + ••• = / — 1. The inverse coefficients 
ai can be obtained by doing a commutation 6/ <H» a/. 
Using the definition a = j/j CJ we can rewrite Eq. (|10[) 

as 



j(*)=iJl + X> [Wr(t)] 1 



(12) 



1=1 



Eq. p^l) represents the time dependence of the current 
density in a flux relaxation process. Putting t = in 
Eq. ([12p , we obtain an initial current density 



3. Time dependence of the current density in a flux 
relaxation process 

To obtain the time dependence of the current density 
jit), we still need the time dependence of the activation 
energy U r {t). In the flux relaxation process, U r (t) is an 
increasing function of the time t. An early studyi^ has 
shown that, with logarithmic accuracy, U r (t) can be writ- 
ten as 



U r (t) = kTln 1 + 



ti+t 



(6) 



where k is the Boltzmann constant, T is temperature, r 
is a short time scale parameter, and 



U=r[e 



(e u *' kT - l) 



(7) 



is a virtual time interval, during which the activation 
energy increase from to the initial value Ui. 
Combining Eq.© and Eq.©, we have 



~(t) = j2 a i (o-i)', 



i=i 



where 



w r {t) = -kTln 1 + 



U+t 



Inverting Eq.©, we have 

oo 



(8) 



(9) 



(10) 



Jr = Jc 



i+5> 



i=i 



- -kTln ( 1 + - 



(13) 



Inverting Ea.(fI5]l (or substituting Eq.© into Eq.Q), 
we obtain the time parameter ti expressed in terms of 
the initial current density ji 



,-a,(ai-l)VfeT _ 



1 



(14) 



where a, = ji/j c . 

Eq.([2]) and Eq. (fl"2"|) are the main formulas of the new 
polynomial model of the flux relaxation process. The 
fitting of Eq. (fl"2"|) to the experimental data from a ring 
shape BiiSr-iCaCu20%j rX single crystal^ is shown in 
Fig.l. 



3. Unified polynomial model of flux relaxation phenomenon 

We have constructed the polynomial model of the flux 
relaxation phenomenon at the critical current density j c 
in Eq.J2]) and Eq. (fl2"|) . We have also constructed the 
polynomial model of the flux relaxation phenomenon at 
a vanishing current density in Ref. [l] and [H. Here, we 
show that these polynomial models can be unified into 
one theory that can be applied to an arbitrary center 
point. 

Let j s G [0,j c ] be an arbitrary center point. Referring 
to Eq.([T]), we can expand the activation energy of a flux 
relaxation process about the point j s as 



i=i 



U r (j) = U r (j s ) + ^2a l (j-j s ) 1 , 
1=1 



(15) 
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the polynomial model of the flux relaxation phenomenon 
constructed at the critical point (see Eq.® and Eq. ([T2")) ). 



B. Polynomial model of vortex penetration 
phenomenon 

In the vortex penetration process, it is convenient to 
use internal field (penetrated field) instead of current 
density^^ii^ In this section, therefore, we expanded the 
activation energy of the vortex penetration process as an 
infinite series of the internal field. 



FIG. 1: (Color online) Polynomial model of flux relax- 
ation. The scattering points are the time dependence of 
the persistent current induced at 28 K in a ring shape 
Bi2Sr2CaCv>20g+ x crystal. [Ref.ll5l| The solid black line is the 
theoretical fit with Eq.(fl2)l: j(t) = j c {l + ]Tf =1 bi[w r (t)] 1 }. 
The fitting results are: j c = (9.36 ± 0.03) x 10 7 , U = 
0.79 ± 0.04, t = 1.18 ± 0.01, 6i = (0.2963 ± 0.0005)/(28fc), 
b 2 = (0.0341 ±0.0001)/(28fc) 2 , b 3 = (0.0014±0.0000)/(28fc) 3 . 



where U r (j s ) is the activation energy at the center point 
j s and ai = Ur\j s )/l\. Here, we cannot use the variable 
j/j s (see Eq.©) because it cannot be applied to the zero 
point j s =0. 

Combining Eq.© and Eq. (fl"5|) . we have 



i=i 



where 



w r (t) = kTln 1 + 



U + t 



Urijs). 



(16) 



(17) 



The relation between ai and bi is the same as that of a; 
and b[ (see Eq. tfTTj) ). Using Eq.(fT6"|). one can analysis the 
experimental data by choosing an arbitrary center point 

Js 6 [0,Je]. 

Similar to Eq. (|13l) . we have the initial current density 



i=i 



kTln [ 1 + - ] - Urijs) 



(18) 



1. Internal field dependence of the activation energy in a 
vortex penetration process 

In the vortex penetration process, under an applied 
field B a , the internal field B approaches a maximum 
value B e (equilibrium internal field) when the super- 
conductor approaches an equilibrium state. Due to the 
surface screening effect (or Meissner effect), B e must be 
smaller than B a , that is, B e < B a . (Ref. y) 

Furthermore, in the vortex penetration process, the 
activation energy U p is an increasing function of the in- 
ternal field B. To construct a polynomial model about 
the equilibrium internal field B e , let us define a normal- 
ized internal field (3 — B/B e . We have (3 £ [0, 1] because 
Be[0,B e ]. 

Referring to Eq.([T]) and rewriting the activation energy 
as U p [l + i/3 — 1)], we see that U p needs to be expanded 
as an increasing function of (/3 — 1) about the point /3 = 1 
(B = B e ). Let us explicitly write out the series expres- 
sion of U v as 



U p (p) = U e + Y i c l (p-l) 1 



(20) 



i=i 



and the virtual time interval 



where U e = U p (l), c x = U p {l), c 2 = E£'(l)/2!, 

Q — Up il)/ 11. The parameter U e is the equilibrium 
activation energy (at the equilibrium internal field B e ) 
and is also the maximum activation energy in the vortex 
penetration process under the applied field B a . 

Under a vanishing internal field, a vortex inside the su- 
perconductor is subjected to the Bean-Livingstor>ii sur- 
face barrier Ubl and pinning potential U c . Therefore, the 
vortex activation energy at the vanishing internal field 
can be written as 



U 



t < e 



U r (js)/kT 



n 

i=i 



(19) 



The current density j s in Eq.(fT5|) and Eq. fTTjf is an 
arbitrary center point in the interval [0, j c ]. If we choose 
j s = 0, then we obtain the polynomial model of the flux 
relaxation phenomenon constructed at the zero point (see 
Ref. 12 and 3). If we choose j s = j c , then we obtain 



Upo = U BL + U c 



(21) 



On the other hand, U p $ can also be obtained by putting 
P = in Eq.®, that is, 



U 



p0 



Cl, 



(22) 
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Combing Eq.(|2T)l and Eq.(|22]). we have 

oo 

U c = U e + J2(-l) l ci~U BL . (23) 



z=i 



The cocfBcients U e and q can be obtained from exper- 
iments (see Eq. ([28|) ). Doing a measurement at a position 
away from the surface where Ubl = 0, one can calculate 
the pinning potential using Eq.(| 



2. Time dependence of the internal field in a vortex 
penetration process 



To obtain the time dependence of the internal field 
B(t), we still need the time dependence of the activation 
energy U p {t). An early study^ has shown that U p {t) can 
be written as 

U p (t) = kTln (l + \ = kTln (l + to +^ +t 

(24) 

where r is a short time scale parameter and 

U = r(e^/ fcT -l), (25a) 



to = r (e u *°' hT _ t 
t v = ti - t a = t (e 



Ui/kT _ p U p0 /kT 



(25b) 
(25c) 



where to is a time parameter equivalent to the potential 
U p o and t v is a virtual time parameter, during which the 
activation energy increase from U p q to the initial value 
Ui. 



Combining Eq.(|20|) and Eq.pljl. we have 

oo 

UJp(t) = 5^ciC8-l)', (26) 
i=i 

where 

w p (t) = kTln (\ + t o + t v+t \ _ Ue (2?) 



Inverting Eq. ([2^1) and using the definition j3 — B/B e , 
we have 



B(t) = B e h+J2 d i KW] 



(28) 



1=1 



The coefficients c; and di can be obtained by doing a 
replacement ai — > q and bi di in Eg. (1111) . 



Eg . ([28)1 represents the time dependence of the internal 
field in a vortex penetration process. Putting t = in 



Eq. (|28]) . we obtain an initial internal field 



B i = B e ll + J2d l 



i=i 



1 1 



(29) 

If the initial internal field B l = 0, then t v = 0. Eq. (|2^|) 

becomes 



kTln 1 



= 1 + d - 
i=i 

oo 

= 1 + di [U pQ - U e ] 



tr 



u, 



(30) 



i=i 



Eg . ([3^| can also be obtained by inverting Eq. 
Substituting Eq.® into Eo. (j2"5aj) , we have 



/, = r^e tr '/ fcr n eC,(A ~ 1)VfcT - 1 



(31) 



where ft = Bi/B e . 

Substituting Ea.(|20]l into Eq. (|25cp . we have 



t v = t S e L 



u./trrr qft-ij'/tT _ p u p0 /kT 



n 

Z = l 



(32) 



Eq. ([20)) and Eq. ([28|) are the main formulas of the new 
polynomial model of the vortex penetration process. The 
fitting of Eq. (|2"51) to the experimental data from a disk 
shape BiiSriCaCu20%+ x single crystal^ is shown in 
Fig.2. 



3. Unified polynomial model of vortex penetration 
phenomenon 



We have constructed the polynomial model of vortex 
penetration phenomenon at the equilibrium internal field 
B e in Eq. ([2U| and Eq. ([2"5|) . We have also constructed the 
polynomial model of vortex penetration phenomenon at 
a vanishing internal field in Ref. |2j and |3|. Similar to 
Eq. ([T5)) and Eq.fllT)]). let us now construct the unified 
polynomial model of vortex penetration phenomenon at 
an arbitrary center point. 

Let B s G [0, B e ] be an arbitrary center point. Referring 
to Eq. ([T]) , we can expand the activation energy of a vortex 
penetration process about the point B s as 



U P (B) = U P (B S ) + J2 HB - B s ) 1 



(33) 



i=i 



where U P (B S ) is the activation energy at the center point 
B s and ci=U^{B s )/l\. 
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FIG. 2: (Color online) Polynomial model of vortex pen- 
etration. The scattering points are the time dependence 
of the internal field induced at 28 K in a disk shape 
Bi2Sr2CaCv>20g+ x crystal. [Ref.ll5l| The solid black line is the 
theoretical fit with Eq.((28j): B(t) = B e {l + £? = i di[w p (t)] 1 }. 
The maximum internal field B e = (1.36 ± 0.01) x 10" 4 T is 
determined by the experimental measurements. The fitting 
results are: r = 0.019±0.005, t v = 0, to = 0.006±0.001, U e = 
(12.641 ± 0.001) x (28fc), d\ = (0.022 ± 0.000) x 10" 4 /(28fc), 
d 2 = -(0.0049 ± 0.0008) x 10~ 4 /(28£;) 2 , d 3 = (0.00125 ± 
0.00006)/(28fe) 3 . (The B e should be multiplied by a geome- 
try factor because the field B(t) is measured with a magnetic 
sensor positioned above the sample. But this does not prevent 
us from finding out other parameters.) 



In Eq.(ja3J) and Eq.([34]), if we put B s = 0, then we 
obtain the polynomial model of the vortex penetration 
phenomenon constructed at the vanishing internal field 
(see Ref. [U and [H). If we put B s = B e , then we ob- 
tain polynomial model of the vortex penetration phe- 
nomenon constructed at the critical point (see Eq.(|20p 
and Eq.(g5|). 

Finally, it should be mentioned that the "center 
points" in the unified polynomial models are not ob- 
tained by data fitting, but were chosen by the authors 
in advance. In other words, the center points are given 
numbers. This means that, in a vortex penetration pro- 
cess, one needs to know either the equilibrium internal 
field B e or equilibrium activation energy U e in search- 
ing for the physical parameters at the critical point (see 
Fig. 2). B e is determined by the applied field B a , i.e., 
B e < B a . In a flux relaxation process, however, the ac- 
tivation energy U r is zero at the critical current density 
j c . Therefore, one can directly apply Eq. (fT2"|) or Eq. fTfi 
(put j s = j c ) to experimental data to obtain j c . 



IV. INVERSE MODEL 

In this section, we first presented a generalization to 
the inverse models of the flux relaxation phenomena dis- 
cussed in Sec. HTJ Next, we applied this method to a 
vortex penetration phenomenon. 



Combining Eq.([M| and Eq.([33]), we have 

oo 

B(t)=B s +J2dl[wp{t)} 1 > (34) 

where 



i=i 



W p (t) = kTln 1 + 



- U P (B S ). (35) 



A. Inverse model of flux relaxation phenomenon 

In the inverse model of a flux relaxation process, the 
activation energy is supposed to be expanded as an in- 
finite series of the inverse current density 1/j about the 
critical current density j c . 



Current dependence of the activation energy in a flux 
relaxation process 



The relation between q and di is the same as that of 
ai and bi (see Eq. (fTT|) ). 

Similar to Eq. (|2T))) . we have the initial internal field 

-i l 



Bi = B S +J2 di 



i=i 



kTln 1 



t, 



U P {B S ) 



the virtual time interval 



U p (B.)/kT JJ e ci(Bi-B,) l /kT 

i=i 



and 

{oo 
e U p (B e )/kTTT e ci(Bi-B 3 ) l /kT _ e U p0 /kT 
1=1 



(36) 



(37) 



(38) 



In the flux relaxation process, the activation energy U r 
is a decreasing function of the current density j. There- 
fore, U r must be an increasing function of To con- 
struct the inverse model, let us define a new variable 
^ = 3c/ j, where j c is the critical current density. We 
have A G [1, oo) because j G [0,j c ] (see Sec. El- 
Referring to Eq.([T]) and rewriting the activation energy 
as U r [l + (A — 1)], we see that U r needs to be expanded 
as an increasing function of (A — 1) about the point A = 1 
(j = jc)- Let us explicitly write out the series expression 
of U r as 



U r (X) = J2 ei {X-l) 1 



(39) 



i=i 



where U r (l) 
Ul'(l)/2\,.. 



= (U r is zero at j c ), e\ 
e l = uP(l)/l\. 



K(l), e 2 
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At a first look, one may doubt the correctness of 
Ea.([55]l. To remove this doubt, let us check Eq.® with 
several existing activation energies. 

a. Inverse-power activation energy®' 7 



U 



p 1=1 
b. Logarithmic activation energy&& 



(ji - 1) ■ ■ ■ Qti - I + 1) A ; 



U ln\ = U Y^ 



1 = 1 



c. Linear activation energy 11 



Eq. (l4Tj) . we obtain an initial current density 



Jt =Jc 



i=i 



kTln I 1 + - 



(42) 



Inverting Ea.(|l2"]) (or substituting Eq.® into Eq.©), 
we have 



u=r m e ei(A *- i)!/fcT -i , 



(43) 



where A l = j c /j;. 

Eq. (|39|) and Eq. ([4"Tj) are the main formulas of the in- 
verse model of the flux relaxation process. 



f»(i-{)=^i;tiP(A-i)' 



1=1 



The other activation energies^ can also be checked in 
the same way, as desired. This confirms that the ex- 
isting inverse activation energies are the special cases of 
Eq. (f39f . Now it is clear that the activation energy U r can 
be expressed as an infinite series of the inverse current 
density 1/j. 

Finally, it should be mentioned that the pinning po- 
tential U c needs to be discussed at a vanishing current 
density (J = 0). But in the inverse model, j = is a 
singularity of U r (see Eq.(J35])). Thus, we do not discuss 
U r in the inverse model. 



2. Time dependence of the current density in a flux 
relaxation process 

Combining Eq.© and Eq . (|39|) . we have 

oo 

-uv(i)=5>(A-l)<, (40) 
i=i 

where w r (t) is defined by Eq.©. 

Inverting Eq. (j40| and using the definition A = j c /j, we 
have 



j(t)=jc\l+J^fl[-W r (t)} 



(41) 



i=i 



The coefficients e; and /; can be obtained by doing a 
replacement ai — > e; and bi — > /; in Eq. ([TT]) . 

Eq. (|41[) represents the time dependence of the current 
density in a flux relaxation process. Putting t = in 



Inverse model of vortex penetration 
phenomenon 



In the inverse model of a vortex penetration process, 
the activation energy is supposed to be expanded as an 
infinite series of the inverse internal field 1/B about the 
equilibrium internal field B e . 



1. Internal field dependence of the activation energy in a 
vortex penetration process 

In the vortex penetration process, the activation en- 
ergy Up is an increasing function of the internal field B. 
Thus, Up must be a decreasing function of 1/B. To con- 
struct the inverse model, let us define a new variable 
fj = B e /B, where B e is the equilibrium internal field. 
We have a G [1, oo) because B e [0, B e ] (see Sec. N\i . 

Referring to Eq. |T]) and rewriting the activation energy 
as Up[l+(a— 1)], we see that U p needs to be expanded as a 
decreasing function of (a— 1) about the point a = 1 (B = 
B e ). Let us explicitly write out the series expression of 
U v as 



U„ 



u e 



9l{<r-l) 1 



(44) 



where U e = U p (l), g x = -U' p {V), 



91 



U^ l \l)/l\. The[/ e inEq. 



.g 2 = -rj;'(i)/2!, 

is the equilibrium ac- 



tivation energy corresponding to the equilibrium internal 
field B e (cr — 1). This is already discussed in Eq. (|2"Ul) . 

The pinning potential U c needs to be discussed at a 
vanishing internal field (B = 0). But B = is a singu- 
larity of Up (see Eq. (|4"4")0 . Therefore, we do not discuss 
U c in the inverse model of the vortex penetration process. 



8 



2. Time dependence of the internal field in a vortex 
penetration process 

Combining Eq. ([M|) and Eq.([4"I)). we have 

oo 

-w p {t)=Y J 9l{v-l)\ (45) 
i=i 

where w p (t) is defined in Eq.tPZT]). 

Inverting Eq. p5)l and using the definition a — B e /B, 
we have 



(46) 



i=i 



The coefficients gi and hi can be obtained by doing a 
replacement ai — > gi and bi — > hi in Eq. (|lll) . 

Eq. (|46[) represents the time dependence of the internal 
field in the vortex penetration process. Putting t = in 
Eq. ([46[) . we obtain an initial internal field 



B, = B J 1 + J2 hi 



i=i 



U e - kTln 1 



U 



(47) 



Inverting Eq.fllTl) (or substituting Eq.(g3]) into 
Eq. (|25ap ). we have 




where ai = B e /Bi. 

Eq. (|4~4")) and Eq. (l4T>l) are the main formulas of the in- 
verse model of the vortex penetration process. 



DISCUSSION 



according to the physical parameters that he/she wishes 
to find out from the experiments. 

b. Comparison of polynomial models and inverse 
models In the polynomial models, the activation en- 
ergies are well defined functions in the entire domain 
(Eq.@ and Eq. (p?0)) ). Using these formulas, one can cal- 
culate various physical parameters in the flux relaxation 
process and vortex penetration process, such as the pin- 
ning potential, critical current density and equilibrium 
internal field. 

In the inverse models, however, the activation energies 
diverge at the zero points. This defect limits the appli- 
cation of the inverse models. Let us now discuss them 
separately. 

In the inverse model of a flux relaxation process, the 
activation energy Eq. (|39| diverges at j = 0, although 
it is a generalization of other inverse current dependent 
activation energies. Thus, the inverse model cannot de- 
scribe the vortex behavior at j = and is only valid at 
the region close to the critical current density j c . Also, 
it cannot fit to the experimental data of the entire flux 
relaxation process, because j reduces from j c to with 
increasing time t . This has been confirmed in the fit of ex- 
perimental data (It is not shown here) . Using the inverse 
model, therefore, one meets difficulty in determining the 
physical parameters at the zero point, such as the pinning 
potential U c . Let us now estimate the smallest current 
density that can be described by the inverse model. 

Consider that the maximum activation energy in a 
flux relaxation process is the pinning potential U r Q (see 
Eq.©). Using Eq.®, we have, 



U r (X) =^e,(A-l) ; < U r0 . 
i=i 



Simplifying this inequality, we have 



On the basis of present work and the early worksi~— , 
we can now summarize the theoretical models of vortex 
dynamics in Table [IJ 

a. Comparison of the polynomial models constructed 
at zero points and critical points The physical param- 
eters of a flux relaxation (vortex penetration) process 
can be obtained using a polynomial model constructed 
about either the zero point or the critical point. But 
from the data fitting, we see that the polynomial models 
constructed at the critical points have the advantage in 
determining the physics parameters at the critical points, 
such as the critical current density j c in the flux relax- 
ation process, and the equilibrium internal field B e (or 
equilibrium activation energy U e ) in the vortex pene- 
tration process. However, the polynomial models con- 
structed at the zero points^— have the advantage in de- 
termining the physics parameters at the zero points, such 
as the pinning potential U c . Thus, for the convenience of 
calculation, one should choose a proper theoretical model 



=i 



(49) 



Thus, in the inverse model of a flux relaxation process, 
Eq. (l4"Tj) can only be applied to a system with a current 
density satisfying the inequality in Eq. (|49[) . 

Similarly, in the inverse model of a vortex penetration 
process, the activation energy Eq. (l4"4")l diverges at B = 0. 
Thus, it cannot describe the vortex behavior at B = 
and is only valid at the region close to the equilibrium 
internal field B e . Also, it cannot fit to the experimental 
data of the entire vortex penetration process, because B 
increases from to B e with increasing time t. Therefore, 
one meets difficulty in determining the physical param- 
eters at the zero points using the inverse model. Let 
us now estimate the smallest internal field that can be 
described by the inverse model. 

Consider that the minimum activation energy in a vor- 
tex penetration process is the pinning potential U p o (see 
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TABLE I: Vortex activation energies in flux relaxation and vortex penetration process. 



Flux Relaxation 



Vortex Penetration 



I a Infinite scries [Refs. 1 3] 

Quadratic [Refs. h-S\ 
Linear [Refs. ^,-4] 



U a (j) = U c -Y,a i j i , 
i=i 

U a (j) = U e - axj - a 2 j 2 , 
U a (j) =U C - aij 



U b {B) = U c +Y,aiB l 1 

U b {B) = U c + ai B + a 2 B 2 , 
U b {B) = U c + aiB 



II 6 


Infinite series 


U r (a) = - J2 a i ( a - 1) ) 


U P (J3) = 


U e + 






Quadratic 
Linear 


U r (a) = —ai (a — 1) — a 2 (a — 1) , 
U r (a) = — ai (a — 1) 


U p {[3) = U e 4 

u P W) 


-ciC8 

= c/ e 


- 1 l) + c 2 (/3-l) 2 , 
fci(^-l) 



III c Infinite series 

Inverse-power [Refs. ,6,7] 
Logarithmic [Refs. .8,9] 



U r (X) 



E e;(A - 

(C/ /m)"(V - 1) 
UqIuX 



I) 



U p {a) =U e - E 5/(^-1) 

i=r 



"Polynomial models constructed about zero point. 

'Polynomial models constructed about critical point, a = j/j c , 
= B/B e . Expanding the quadratic and linear activation energy, 
one find that they have the same form as those in I. 

c Inverse models constructed about critical point. A = j c /j, a = 

B e /B. 



Eq.([2l]>). Using Eq.([i4]). we have 



U p {a) = U e -Y J 9i{° ~i) l >Uj 



pO- 



i=i 



Simplifying this inequality, we have 



B(t) > B e 



l + J2hl{U e - U p0 ) 1 



1=1 



(50) 



Thus, in the inverse model of a vortex penetration pro- 
cess, Eq. (|46|) can only be applied to a system with an 
internal field satisfying the inequality in Eq. (|50[) . 

c. Surface barrier The Bean-Livingston surface 
barrier— affects vortex motion. It is usually considered 
along with the pinning potential at the zero points. In 
the polynomial models, we directly combined the sur- 
face barrier into the activation energies (see Eq.([3]) and 
Eg . (l2lT) ) . In the inverse models, however, the activation 
energies diverge at the zero points. Thus, we do not con- 
sider the surface barrier in the inverse models. 

We have ignored the geometrical barrie r 12 ' 13 in the 
derivation, because it depends on the sample shape. This 
effect is significant only for some special geometrical 
shapes. 

d. Anisotropy of superconductor The layered super- 
conductors are strongly anisotropic. This property can 



be represented by the coherence length and penetration 
depth. Because the activation energies are the func- 
tions of the coherence length and penetration depth, the 
anisotropy of the superconductors expresses itself in the 
flux relaxation and vortex penetration process. In this 
sense, we expect the derived equations to be applicable 
to both low-T c and high-T c superconductors, although 
we did not explicitly consider the anisotropy of the su- 
perconductors in our derivation. 



VI. CONCLUSION 

The polynomial models of flux relaxation and vortex 
penetration phenomena can be constructed about the 
critical points. These models are well defined functions in 
the entire domain. They have advantage in determining 
the critical current density and equilibrium internal field, 
comparing with the early proposed polynomial models 
which were constructed at the zero points. The polyno- 
mial models constructed at the zero points and critical 
points can be unified into one theory. On the other hand, 
the inverse model of the flux relaxation and vortex pen- 
etration phenomena were also constructed at the critical 
points. The inverse models diverge at the zero points and 
are only valid at the region close to the critical points. 
They cannot be used to determine the pinning potential. 
Therefore, the polynomial models have advantages and 
give better description to the vortex dynamics. 
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